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ABSTRACT 

Current theories of structure formation predict specific density profiles of galaxy dark matter 
haloes, and with weak gravitational lensing we can probe these profiles on several scales. On 
small scales, higher-order shape distortions known as flexion add significant detail to the weak 
lensing measurements. We present here the first detection of a galaxy-galaxy flexion signal in 
space-based data, obtained using a new Shapelets pipeline introduced here. We combine this 
higher-order lensing signal with shear to constrain the average density profile of the galaxy 
lenses in the Hubble Space Telescope COSMOS survey. We also show that light from nearby 
bright objects can significantly affect flexion measurements. After correcting for the influence 
of lens light, we show that the inclusion of flexion provides tighter constraints on density 
profiles than does shear alone. Finally we find an average density profile consistent with an 
isothermal sphere. 
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1 INTRODUCTION 

Weak gravitational lensing is a powerful technique for studying the 
distribution of matter in the universe due to its ability to model 
the matter distribution in foreground structures, independent of 
the nature of the matter present. As the light from background 
■ sources is bent around foreground lenses, the galaxy images get 
' distorted by the tidal gravitational field. The first-order distortion 
is known as shear and is essentially an elongation of the image 
causing the source galaxy to appear stretched in one direction. 
This type of distortion measurement has been used in a wide vari- 
ety of cosmological studies ranging fro m modeling the large-scale 
struct u re using cosmic shear (see e.g. IVan Waerjjeke & Melliej 
' l2003l : iHoekstra & JainI l2008l : iMunshi et all l2008l . for reviews) 
to determining galaxy h alo shapes u sing galaxy-galaxy lens- 
ing (Jioekstra et al. 20041 : iMandelbaum et al. 2006a: Parker et al.^ 

l2007h . 

First described bv lGoldberg & NataraianI (|2002|) . the second- 
order distortio n is a relatively new addi t ion which has sin ce been 
named flexion ( Goldber g & Baconll2005l : lBacon et al I2OO6). There 
are two types of flexion relevant to weak lensing studies: the 
first flexion induces a skewness of the brightness profile whilst 
the second flexion is a three-pronged distortion. In combination 
with shear, these distortions cause the well-known banana shape of 
lensed source images. As flexion is effectively the gradient of shear, 
it is sensitive on small scales. This makes it an important comple- 
ment to shear which is sensitive on relatively large scales only. By 
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virtue of this, and of the orthogonality of the three measurements, 
flexion is highly beneficial to investigations of the inner profiles of 
dark matter haloes, where baryons become important, and to the 
detecti on of substruc ture in cluster haloes. Indeed, it was recently 
shown ter et al.ll201Ql) that mass reconstructions profit from the use 
of flexions in combination with shear, and flexion has already been 
used to constrain the halo mass distribution and to detect substruc - 
ture in clusters of galaxies jLeonard et al .l201Ql : IOkura et al.l2008l) . 
To provide more information on substructure and mass profiles, 
there are currently new statistical flexion tools being developed (eg . 
[Leonard et aP l2009l : [Leonard & Khi g '2010"; 'Bacon et al. '2010). 
Another application, as discussed in[Hawken & Bridle ( 2009), is to 
use both flexions in combination with shear to significantly tighten 
the constraints on galaxy halo ellipticities compared to using shear 
alone. 

The shape measure ment tech nique known as Shapelets 
(lRefregiejr2003 : Refregier & Bacon|[2003.) works by decomposing 
a galaxy image into a series of 2D Hermite polynomials. These 
provide a simple framework for describing the main galaxy image 
distortion operators, such as shear and flexion, and the convolution 
with the point- spread function (PSF). Due to the flexible treatment 
of the PSF, the Shapelets formalism has an advantage over the cur- 
rently most widely used shape mea surement method, KSB (from 
iKaiser. Squires. & Broadhurs3ll995h . since KSB uses an idealised 
model for the PSF whilst Shapelets is more versatile. The KSB 
equivalent for flexion is known as HOLICs dOkura et al.ll2007h . 

Since the field of weak lensing is relatively new, lensing mea- 
surements are continuously being improved in accuracy and ap- 
plicability. Being a statistical technique, however, the accuracy 
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of the weak lensing results depends heavily on the amount of 
data available. G alaxy-galaxy flexion has been tentatively observed 
(iGoldberg & Baco n 2005) using the ground-based Deep Lens Sur- 
vey (DLS), but to further investigate galaxy-size haloes more and 
better data is needed. With large surveys such as the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) and the Red 
Sequence Cluster Surveys (RCS, RCS2) available, and new sur- 
veys like the 1500 square degree Kilo-Degree Survey (KiDS) im- 
minent, the future looks bright. However, a space-based data set 
provides better resolution and such a data set is already accessi- 
ble to us: the HST COSMOS survey. Using this data we will in 
this paper improve on the galaxy-galaxy flexion measurements of 
iGoldberg & BaconI (l2005h . 

This paper is organised as follows: in Section[2lwe review the 
formalism for shear and flexion, whilst we review the Shapelets 
method in Section [3] with a description of our implementation 
(dubbed the MV pipeline) in Section [TT] In Section |4] we test the 
MV pipeline on simulations and in Section[5]the pipeline is applied 
to data from the COSMOS survey. We conclude in Section [6l 

Throug hout this paper we assume the following cosmology 
fWMAP7: IKomatsu et ahkOl O*): 

{QM,^A,h,a8,w) = (0.27,0.73,0.70,0.81,-1) 



2 SHEAR AND FLEXION 

We begin by briefly reviewing the weak lensing formalism. 
Flexion is a second-ord er lensing effect first int roduced by 
'Goldberg & Bacon (2005) and further developed bv 'Bac on et al.l 
(2006) (hereafter B06). It arises from the fact that convergence and 
shear are not constant across a source image, and can be used to 
describe how these fields fluctuate. In the weak lensing regime, the 
lensed surface brightness of a source galaxy, /(x), is related to the 
unlensed surface brightness, /o(x), via 



/(x) 



1 + 



(A I)ijXj + —DijkXjXk 



(1) 



Here / is the identity matrix, Xi denotes lensed coordinates, and 



At - 71 -72 
—72 1 — At + 71 



(2) 



with 



a second derivative of the lens- 



ing potential ^, where subscripts denote partial differentiation. 

71 = I (V^cccc — i/^yy) and 72 = i/^xy 2ltq thc two componcnts of the 
complex shear 7 = 71 + 272. The matrix 



D 



ijk 



dxk 



(3) 



describes how convergence and shear vary across a source im- 
age. We can re-express Dijk as the sum of two flexions: 
Dijk = J^ijk + Gijk- The two flexions, the first flexion T (known 
as F flexion or one-flexion) and the second flexion Q (known as 
G flexion or three-flexion), are the derivatives of the convergence 
and shear fields. There are four flexion components, each of which 
may be written in terms of the third derivatives of the lensing po- 
tential dHawken & Bridle 2002): 



2 {^'l/^xxy 



(J) 



where = J^i + 1^2 and Q = Qi-\- iQ2 are the complex F and G 
flexions respectively. The full matrices J^ijk and Gnk in terms of 
the four flexion components are written explicitly in B06. 



3 SHAPELETS 



The S hapelets basis function set was introduced by iRefregiej 
(l2003h and is more fully described there. In summary, the surface 
brightness of an object /(x) can be expressed as a sum of orthogo- 
nal 2D functions 



(8) 



where Sab are the Shapelets coefficients and the Shapelets basis 
functions Bafe(x; /3) are defined as 



Bab{^;l3) = kabr''e ^ Ha{x/ 0)Hb{y/ 0). 



(9) 



Here kab is a normalization constant, [3 is the Shapelets scale ra- 
dius, {x,y) are coordinates on the image plane and Hn{x) is a 
Hermite polynomial of order n. The Shapelets basis functions are 
easily recognised as the energy eigenstates of the 2D Quantum Har- 
monic Oscillator (QHO). The formalism developed for the QHO 
can also be applied to Shapelets, providing analytical expressions 
for transformations such as shear and flexion. In theory, an object 
can be perfectly described through a decomposition into Shapelets 
up to order n ^ 00 but in practice the expansion has to be trun- 
cated. We truncate at combined order rimax = a + 6 to avoid intro- 
ducing a preferred direction. 

Convolution with the point- spread function (PSF) can also be 
done analytically in the Shapelets formalism by simply multiplying 
the Shapelets expansion by a PSF matrix P: 



a\a2h\b2 



(/^obj , /^c, 



'/3con/3obj(^psf ^/5con/3objy^psf (^(\\ 



where pab are the Shapelets coefficients of the PSF and /3psf , /3obj 
and [3 con are the scale radii of the PSF, the object and the resulting 



PSF convolved object respectively. C!^^i^ is a convolution tensor 
which d epends on the dif ferent scale radii and the full expression is 
given in lRefregied (l2003h . 



3.1 The MV Pipeline 

We introduce here an implementation of the Shapelets method 
which builds on a previ ous implementation described in iKuiikenI 
(l2006h rhereafter lKK06h . This approach creates a Shapelets repre- 
sentation of the brightness profile of a PSF-convolved galaxy im- 
age. It also creates a model circular source and applies shear and 
flexion to it before convolving it with the point- spread function 
(PSF). Finally it fits the galaxy image to this modeled source in 
order to find the amount by which it has been sheared and flexed. 

To first order in ellipticity s and flexions / and g, the model 
object can be written as 



Gi 



2 (V^cccccc H~ "^^yyx^ 
— {'ipxxy + ^^yyy) 



XXX 



(4) 
(5) 
(6) 



^c„C" (11) 



where P is the PSF matrix, T\ S\ and G' are the transla- 
tion, shear, F flexion and G flexion operators respectively and ti,Si, 
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Figure 1. Polar Shapelets basis functions up to a maximum Shapelets order 
of rimax = 10. For m ^ 0, the real components of the basis functions are 
shown whilst for m < the imaginary components are shown. The solid 
purple (thick) lines mark the coefficients used by the MV pipeline to esti- 
mate the shear and flexions for an analysis with rimax = 10. The dashed 
purple (thick) lines mark the coefficients not used by the iKK06 implemen- 
tation for the same rimax- 



fi and gi are the corresponding coefficients. The translation terms 
here ensure that fits spoiled by undue centroid shifts are caught. The 
operators are acting on a circular source which can be expressed as 
a series of circular Shapelets with coefficients Cn where n is 
even and the series is truncated at Nc — rimax — 2. The reason for 
truncating at Nc rather than rimax is to safeguard against PSF struc- 
ture at higher orders affecting the highest order Shapelets used. To 
avoid introducing signal-to-noise (S/N) dependent biases, the rimax 
is kept constant for all galaxies rather than being allowed to vary 
according to size or brightness. For faint sources, this means the 
higher-order coefficients will be noisy but unbiased. 

Once we have a cartesian Shapelets representation of both the 
sheared, flexed and PSF convolved circular model and of the PSF 
convolved object we want to fit, we convert them both into po- 
lar Shapelets as described in iKKOd Polar Shapelets are simply 
cartesian Shapelets of order n = a + 6 expressed in polar coor- 
dinates, resulting in polar Shapelets of order n with angular order 
n and n + m even. The construction of these is discussed 
in lRefregiei](l2003h and further in vestigated in Massey & Refregier 
(l2005h and Massev et al.l ( [2007bh . In our implementation, the pur- 
pose of converting the model and object Shapelets expansions into 
polar Shapelets is to avoid truncation effects. F flexion, shear and 
G flexion operators acting on a polar Shapelet of order (n, m) gen- 
erate terms at order (nil, m=bl), (n±2, m±2) and (n=b3, m±3) 
respectively. By truncating the polar Shapelets expansion in the di- 
amond shape shown in Figure [T] i.e. only including terms up to 
order {Nc, 0), {Nc - 1, ±1), (A^e - 2, ±2) and (A^e - 3, ±3) in 
the fit, we minimise truncation effects from the mixing of orders. 

The model is fit to each source using least- squares, resulting 
in a simultaneous estimate for the ellipticity (si, the F flexion 
(/i 5/2), and the G flexion {gi , ^2)- As explained in "KK06", the er- 
rors on the Shapelet coefficients are derived from the photon noise 
and propagated through the function for this fit. By differentiat- 



ing the at the best-fit, we obtain the covariances between the fit 
parameters, resulting in proper error estimates. 

In essence, the main development since KK06 is the addi- 
tion of flexion to the model and the inclusion of higher order polar 
Shapelets (m = ±3) in the fit. 



4 TESTING THE PIPELINE 

Several aspects of the pipeline, such as the choice of scale radius /3, 
the method of PSF correction and the effect of noise on ellipticity 
estimates, have b een tho roughly tested in KK06 as part of the de- 
velopment of the lKKOd pipeline. In this section we will therefore 
focus on testing the recovery of shear and flexion. 



4.1 GREAT08 

As participants in the GRav itational lEnsing Accura cy Testing 
2008 (GREAT08) challenge ( Bridle et all [200^ . l201Qh . w e were 
able to contrast the shear measurement capability of the iKKOd 
pipeline with that of the MV pipeline under different observing 
conditions. The challenge provided a large number of simulated 
sheared and pixelated galaxy images with added noise. The per- 
formance of the different shape measurement pipelines taking part 
was quoted in terms of a quality factor, or Q-value, defined as 



Q 



(12) 



where cr^ = aJtat + erf 



Stat ^ ^syst is a Combination of the statistical spread 
in the simulations and the expected systematic errors. The super- 
scripts m and t denote measured and true (input) values respec- 
tively and ^ij is the shear component i for simulation image j. 
The differences between the measured and true shears are aver- 
aged over different input shear sets k and simulation branches 
/. The whole expression is normalised by kq so that a method 
with a purely statistical spread in the measured shears will have 
a Q-value of kq which is the level desirable for future surveys. 
In the case of GREAT08, kq = 1000 and = 10"'^, giving 
a Q-value nominator of 10~^. With this definition current meth- 
ods, like those t hat took part in the earlier Shear TEstin g Pro- 
gramme (STEP) (iHevmans et aLll2006l : iMassev et alJl2007ah . gen- 
erally achieve 10 < Q < 100. This is sufficient for current weak 
lensing surveys. For a more in-depth discussion on the Q-value and 
its relation to the STEP p arameters m (multiplicative bias) and c 
(additive bias), we refer to Kitching et al.l (l2008l). 

The overall Q-value was similar for the iKKOd and the MV 
pipelines, both in the LowNoise_Blind competition (Q ~ 20) and 
in the RealNoise_Blind (Q ^ 25). When broken down into the sep- 
arate observing condition branches some differences became ap- 
parent. In general the MV pipeline did exceptionally well under 
"good" observing conditions, e.g. for the high S/N branch or for 
well resolved galaxies. Our own simulations described in the next 
section will further test the dependence of the MV performance on 
different observing conditions. 



4.2 FLASHES 

As there is no flexion simulation set publicly available to date, 
we create our own FLexion And SHEar Simulations (FLASHES). 
FLASHES are very similar to the GREAT08 simulations in several 
respects. First, each galaxy is generated on a grid, ensuring that 
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Table 1. The different branches of FLASHES. Four parameters are varied 
between the branches according to this table. 
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there is no overlap of objects, thus avoiding deblending issues. Sec- 
ond, each simulation image consists of 10000 such objects. Third, 
each galaxy is generated through the following sequence: (i) simu- 
late a sheared and/or flexed (elliptical) galaxy model (depending on 
simulation branch); (ii) convolve with the PSF; (iii) apply the noise 
model. Four parameters are varied between the different FLASHES 
branches; the intrinsic galaxy shape, the light profile of the galaxies, 
the S/N of the galaxies and the shape of the PSF. These parameters 
are detailed below and summarised in Table [T] 



4. 2. 1 Simulation Details 

All parameters except for the intrinsic ellipticities are kept con- 
stant in each simulation image, and all images are created using 
Monte-Carlo s election. This is very similar to the process described 
in|KK06 and in lBridle et al.l fcolO), but with the photon trajectories 
being influenced by flexion as well as by shear if required. 

The galaxies are modeled with Sersic intensity profiles 
(ISersidll968 h with varying indices n. A Sersic in- 
dex of n = 0.5 is a Gaussian profile whilst n — 1 and n = 4 
are exponential and de Vaucouleur profiles respectively. Half of 
the FLASHES branches have intrinsically round galaxies whilst 
the other half consists of galaxies with intrinsic ellipticities picked 
randomly from the ellipticity distribution of objects in the COS- 
MOS survey. There is no intrinsic flexion included. The PSFs ap- 
plied to the simulations are nearly Gaussian with a Moffat profile 
/psF = (1 + r^/a^)""^ of index m = 9. In half of the branches, 
the PSF is round whilst in the other half it is elliptical in the hor- 
izontal direction with ei,PSF = 0.02. To mimic the properties of 
the COSMOS survey, we use a PSF FWHM of 2. 1 pixels and a PSF 
convolved galaxy size of 5.8 pixels which is the typical size of the 
galaxies we use in our COSMOS analysis. Finally there are four 
S/N branches, with S/N being defined as Flux/(Flux error). It is ex- 
pected that shape measurements will be less accurate at low S/N. 
For this reason the MV pipeline applies a S/N cut at 10 in general. 
The low S/N branch of 8 is designed to test how biased measure- 
ments are below this cut. The high S/N branch of 100 tests biases 
under near-perfect noise conditions. 

The strength of the different distortions is picked randomly 
but with the following maximum values: |7i,2| ^ 0.05, |^i,2| ^ 
0.008 pixel" ^ and \Qi,2\ ^ 0.02 pixel" ^. The value of each dis- 
tortion component is kept constant across each image, but differs 
between the 30 images in each set, and between different sets. 



4.2.2 Simulation Results 

To estimate the average distortion on each image we use two differ- 
ent techniques: a weighted average with weights inversely propor- 
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Figure 2. The multiplicative bias on the first component for each of shear, 
F flexion and G flexion. The purple stars represent shear, pink circles rep- 
resent F flexion and green triangles represent G flexion. The symbols and 
solid lines show the weighted averages whilst the dashed lines show the 
CHP average. This is from running the MV pipeline on FLASHES, with 
^max = 10. For the results for m2, ci and C2 please refer to Appendix IaI 



tional to the measurement errors, and Convex Hull Peeling (CHP). 
CHP is an efficient way of eliminating outliers and is essentially 
a 2D median. A convex hull, in the context of a point cloud in 
e.g. the 71 , 72 plane, is the minimal convex set of points containing 
that point cloud. Thus if all the points in this convex set were con- 
nected, a polygon containing the entire point cloud would be pro- 
duced. By peeling away convex hulls, outliers are removed from 
the point cloud and the remaining points may be averaged over to 
produce a mean unaffected by extreme results. This is the averag- 
ing technique we used in GREAT08 where we peeled away 50% of 
the measurements before averaging. 

We employ the parameters m and c as used in STEP 
jHevmans et al.ll20o3 : iMassev et al.ll2007ab to quantify the perfor- 
mance of the software: 



/ measured \ 



-7^ 



input I 

m^7- ^ + d 



(13) 



and similarly for the flexions, where i = 1, 2 represents the shear 
component. A negative multiplicative bias nn thus indicates that 
the distortion is generally underestimated. A systematic offset Ci 
may be caused by e.g. insufficient PSF correction. 

In Figure |2] we show the multiplicative bias of the first compo- 
nent for each of shear, F flexion and G flexion as a function of the 
different simulation branches (please refer to Appendix [A] for the 
remaining bias components). For these results we use a Shapelets 
order of rimax = 10. We use SExtractor (iBertin & ArnoutsI 
to detect the objects in each simulation, which we then split 
into clean star and galaxy catalogues by matching to the input cata- 
logue. We keep all properties apart from the one under investigation 
fixed at a fiducial value to allow for a fair comparison. The fiducial 
simulations in Figure |2]have intrinsically round, high S/N galaxies 
with Gaussian light profiles and a circular PSF. 

From the above figure it is clear that both flexions are likely to 
be underestimated, especially for higher Sersic indices. The bias is 
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Figure 3. The shear (top panel), F flexion (middle panel) and G flexion 
(bottom panel) results from galaxy-galaxy lensing simulations, with and 
without Bright Object Removal (BOR). The black stars represent the tan- 
gential signal without BOR and the green triangles represent the same mea- 
surement corrected using BOR. The dashed pink line is the input signal and 
the purple circles are the cross-signal, which is expected to vanish, for the 
uncorrected measurements. Here, the FWHM of the lens is 14 pixels. Note 
the slight underestimation of the shear, the slight overestimation of the G 
flexion and the massive overestimation of the F flexion in the innermost 
bins when BOR is not applied. 



also strongly S/N dependent, particularly for the F flexion. Thus a 
S/N cut is essential to improve the performance of the MV pipeline, 
but a bias correction should also be implemented. Investigating the 
dependence of m on S/N further, we are able to fit the following 
power-law to our FLASHES results: 



mi, 2 



-a(S/N)- 



(14) 



where a and h are constants as follows: for shear 
(a^,6^) = (6.48, 1.78); for F flexion (a^,6^) = (2.30,0.48); 
for G flexion (ag^bg) = (0.36,0.13). We will apply this bias 
correction to our shape measurements in COSMOS, but since 
FLASHES have been tailored for this particular data the biases 
should be explored further before being applied to other surveys. 



4.3 Galaxy-galaxy Simulations and Bright Object Removal 

At the core of weak galaxy-galaxy lensing is the averaging of 
the signal in rings centered on lenses consisting of single galax- 
ies rather than a galaxy cluster. This type of analysis is robust as 
numerous systematics, induced by e.g. the PSF, cancel out. Differ- 
ent systematics may however be introduced, such as the light from 
the central, often brig ht, lens causin g biases in the shape measure- 
ments as discussed in iRowd (l2008h . To study this possible effect, 
we created simple simulations with sources placed in evenly spaced 
rings around a central lens. Apart from source numbers and posi- 
tions, the simulations were created in the same way as FLASHES. 
The S/N of the images was set to 200 to ensure minimum bias, and 
for the same reason the source galaxies had Gaussian light profiles. 
The size and profile parameters of the lens were varied between 
images. 

The results for a lens with an exponential profile are shown in 



Figure [3] (black stars), where we have used rimax = 10. We recover 
a near-perfect average signal in each source circle far from the lens. 
However, close to the lens the shear and G flexion are slightly af- 
fected, but, more strikingly, the F flexion is severely overestimated. 
The conclusion we draw from this is that bright objects can add 
significantly to the F fiexion signal, due to light 'leaking' into the 
Shapelets fitting radius. This causes the pipeline to detect a source 
light profile that is skewed towards the lens, and interpreting it as 
extra F flexion. 

Our solution is to remove any bright objects sufficiently close 
to the source being fit using a technique we introduce here as Bright 
Object Removal (BOR). Before decomposing a galaxy image into 
Shapelets, we identify any bright objects that could conceivably in- 
trude using selection criteria based on distance between the two 
objects, Shapelets fitting radius of the source, and size and bright- 
ness of the intruding object. We then create Sersic models of the in- 
truding objects using GALFIT (Peng et al. 2002) and subtract these 
models from the Shapelets stamp before doing the fitting. It works 
well in these simulations, provided one is careful with the parame- 
ters given to GALF I T as input. The sky background value given to 
GALFIT is particularly important as a small error in this estimate 
results in postage stamp artifacts when the stamps are subtracted 
from the original image. In Figure [3] we also show the results if 
BOR is switched on whilst the rest of the analysis is kept iden- 
tical to the previous run (green triangles). There is still some ex- 
cess F flexion signal around 44 pixels, indicating that there may be 
some residual light remaining, but this excess is smaller than for 
the uncorrected measurements. This provides a confirmation that 
the measured signal reproduces the input signal well if BOR is ap- 
plied, and no new artifacts are introduced. We note, however, that 
the leaking light does not affect the cross component of the mea- 
surements, with the consequence that this effect cannot be detected 
through the usual systematic checks. 



5 COSMOS ANALYSIS 



iGoldberg & BaconI (l2005h made a first detection of galaxy-galaxy 
flexion using the ground-based DLS, proving that flexion can in- 
deed be detected, but ultimately they were hampered by the small 
size of their sample, the lack of redshifts and the extra blurring 
caused by the atmosphere. Therefore we choose the spa ce-based 
Cosmic Evolution Survey f COSMOS. fScoville et all (l2007i,) ) as the 
first real dataset for the MV pipeline. Thanks to the depth of this 
survey we will have access to more than a thousand times as many 
lens-source pairs as Goldberg & Bacon (2005) did. More than half 
of these have photometric redshifts meaning that the division of the 
sample into lenses and sources will be more accurate. The intention 
is to provide independent confirmation that galaxy-galaxy flexion 
has high enough S/N to be detected, and that the software presented 
in this paper is able to do it. We will also look closer to the lens than 
previous analyses and attempt to combine shear and flexion to give 
constraints on galaxy dark matter halo profiles. 



5.1 The COSMOS Data Set 

COSMOS is to date the largest contiguous field imaged by the Hub- 
ble Space Telescope (HST) with a total area of 1.64 deg^. The 579 
tiles were observed in F814W (I-band) by the Advanced Camera 
for Surveys (ACS) between October 2003 and November 2005. 
Each tile consisted of 4 dithered exposures of 507 seconds each 
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(2028 seconds in total) with about 95% of the survey area benefit- 
ing from the full 4 exposures. 

We use the images reduced bv lSchrabback et alJfeOlOh (here- 
after SIO) and also their catalogues for stars and galaxies, detected 
using SExtractor. There are a total of 446 934 galaxies with 
is 14 < 26.7 in the mosaic catalogue, almost half of which have 
COSMOS-30 photometric redshifts from Ilbert et al. (2009). These 
redshifts are magnitude limited and cover the entire COSMOS field 
to a depth of i+ < 25. 

5.2 Data Analysis 

Galaxy-galaxy lensing is less affected by the problems plaguing 
cosmic shear analyses, since most systematic shape distortions in- 
duced by instruments cancel out when azimuthally averaged. Still, 
we have to be careful not to introduce new systematic effects or bi- 
ases, so co rrecting for the PSF and the charge-tran sfer inefficiency 
(CTI) fe.g. lRhodes et al.ll2007l : lMassev et ahlboiOh is important. 

We use all galaxies with redshifts of z < 0.6 as lenses. At 
higher redshifts the light from the lensing galaxies becomes diffi- 
cult to account for due to the small angular separation on the sky, 
as explained further in Appendix |C] Furthermore, imposing a lens 
redshift cut will ensure that the vast majority of sources are truly 
background objects. 

Our source catalogue is comprised of all objects with a shape 
measurement. We clean this catalogue using a series of conditions 
on size and measured shape, detailed in Appendix IBll the most im- 
portant of which is to remove objects with S/N < 10. Roughly 
two-thirds of the remaining sources have individual COSMOS-30 
photometric reds hifts assigned to them. For the remaining third 
(redshift bin 6 in lSlOh we use the estimated redshift distribution 
employed by'Sld to assign mean angular diameter distance ratios 
(Ds/Dis) to each lens-source pair. We are finally left with 216 873 
sources, corresponding to a source density of ^ 37 arcmin"^. For 
the pairs we use, the median lens redshift is ziens = 0.27 and the 
median source redshift is Zsource = 0.98. 

Despite the excellent space-based resolution, we need to 
correct the galaxy shapes for the instrumental PSF. The ACS 
PSF is known to fluctuate both spatially and temporally (e.g. 
[Rhodes et aD l2007l : ISchrabback et all l2007h . a variation mostly 
driven by changes in telescope focus caused for example by the 
breathing of the telescope. We can map the PSF using stars, but, in 
high-galactic latitude ACS fields typically only ^ 10 — 20 stars are 
present. This number is too low for the standard approach of a poly- 
nomial interpolation. Instead, we closely follow the analysis of lSld 
who conducted a principal component analysis (PCA) of the ACS 
PSF variation as measured in dense stellar fields. Details for the 
Shapelets implementation of PCA may be found in Appendix IB 2 1 

A challenge with using CCD detectors in space is that they are 
not protected by the atmosphere. Exposed, they continuously get 
bombarded by radiation, causing deterioration of the chip surface. 
The imperfections created in this way act as charge traps which 
causes inefficiency in the moving of electrons to read-out. This ef - 
fect is known as CTI fe.g. lRhodes et al.ll200llMassev et al.ll201Qh . 
As the electrons get trapped and then released at a later point, 
charge trails following objects are created in the read-out direc- 
tion, effectively causing a spurious shear signal in that direction. 
Our correction for CTI again closely follows . SIO , who derive para- 
metric corrections for the change in polarization for both galaxies 
and stars. For more details on this correction, please refer to Ap- 
pendix |B3] 

Once corrected, the galaxy-galaxy shear and flexion signals 



are weighted according to the geometric lensing efficiency of each 
lens- source pair. In the case of flexion there is an extra scale depen- 
dence of the sig nal. For t he Navarro-Frenk- White (NFW) profile 
(iNavarro, Frenk. & White! Fl 996 ), the strength of the shear signal 
scales as 

7NFW oc — — (15) 

where Di, Dg and Dig are the angular diameter distances to 
the lens, to the source, an d between lens and source respectively 
dWright & Brainerdll2000h . The flexion signals scale as 

-77 n DfDis 

^NFWjt'NFW oc — (10) 

(lB06h . We therefore weight the signals accordingly, scale them to a 
reference lens and source redshift and compute the weighted aver- 
age in 25 logarithmic distance bins (see Appendix IB4l for details). 
We use a reference lens redshift of zi^rei = 0.27 since that is close 
to the effective median redshift of our lenses, and a reference source 
redshift Zs, ref = 0.98. To estimate the errors on each bin and the 
covariances between them, we use 5000 bootstrap resamples of our 
source catalogue. 

5.3 Results 

The results from our galaxy-galaxy lensing analysis of the full 
COSMOS lens and source sample is shown in Figure]?] In the left 
panel we plot the shear results as a function of physical distance 
from the lens. These results agree very well with those from lsiol 
(see Appendix IdI), providing an independent consistency check. To 
this we fit a Singular Isothermal Sphere (SIS) profile and a ten- 
tative NFW profile. Due to the dependenc e on mass a nd redshift 
of the mean concentration parameter (e.g. iDuffv et aP UoOS), the 
NFW profile is only an indication when the spread in lens masses 
and redshifts is as great as it is in the above sample. Splitting the 
sample up into redshift and/or mass bins would increase the confi- 
dence in the fit, but decrease the S/N of the signals significantly. 

The middle and right panels show the F and G flexion results 
respectively, for the same lenses and sources. The profiles plotted 
here are identical to those plotted in the shear panel but translated 
into predictions for T and as opposed to fitted to the flexion data 
directly. The F flexion has a tendency to be overestimated com- 
pared to the predicted profile from the shear, and we investigate 
this discrepancy further in the following sections. We also note that 
we measure a G flexion that is very noisy and consistent with zero. 
This is most likely caused by lack of information in higher m-order 
Shapelets for fainter sources, and we choose to use only shear and 
F flexion in the continued analysis. 

Also shown in pink circles in Figure]?] is the signal if we ap- 
ply the multiplicative S/N- dependent bias correction implied by 
FLASHES. With this correction, the F flexion signal becomes 
slightly higher. This bias correction is only based on one specific 
set of simulations and is thus rather preliminary; this is also in- 
dicated in the increased size of the error bars. Correcting for the 
morphology-dependent bias requires accurate source morphology 
determination. Using the photometric galaxy type estimates from 
lllbert e t al.(2009) as an indicator of morphology we find that < 5% 
of our source sample consists of likely de Vaucouleur objects. This 
type estimate is not accurate enough to implement a morphology 
bias correction, but simply removing the de Vaucouleur candidates 
we identified makes little difference to our results. It is clear, how- 
ever, that an accurate bias calibration of the flexion amplitude, tak- 
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Figure 4. The galaxy-galaxy lensing results for the COSMOS data, using a maximum Shapelets order of rimax = 10. Black solid points represent the 
tangential signal and green circles represent the cross term. The pink circles represent the tangential signal if we apply the multiplicative bias correction 
implied by FLASHES. Note that the SIS and NFW profiles have been fitted to the shear data and then translated into predictions for and Q curves. 
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Figure 5. Comparison between the galaxy-galaxy shear and flexion signals 
with and without Bright Object Removal, showing the non-zero correction 
to the innermost F flexion bin (corresponding to roughly 40 px in Figure [3). 
Black solid points represent the difference between the signals before and 
after correction, with the F flexion in units of kpc~^, whilst green circles 
represent the cross term. 



ing into account both source S/N and brightness profiles, requires 
further investigation. 

5.4 Removing Bright Objects 

We now explore the tendency of the F flexion points to lie above 
the predicted profiles. As shown in Section 14.31 the shape mea- 
sured may be affected by bright objects nearby. We implement 
BOR in our COSMOS analysis to see the effect on real data. For 
very well resolved objects, prominent spiral arms and other com- 
plications cause GALFIT to reject the single Sersic profile fit. Re- 
moving these objects, and the residual light from the wings of the 
profile (Figure [3]), requires a more sophisticated model. For now 
we are only interested in a rough indication of the impact this light 
leakage has on a galaxy-galaxy signal so we will not correct for the 
few large objects in this paper. However, as shown in Figure [S] the 
correction to the innermost F flexion bin is non-zero even without 
accounting for the very large objects. The shear is largely unaf- 
fected, but for flexion analyses in future deeper and larger surveys 
it will be important to correct for this effect. 
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Figure 6. Simulated shear and flexion signals azimuthally averaged in 
galaxy haloes with and without TSIS subhaloes. Grey stars, circles and 
triangles represent the binned shear, F flexion and G flexion respectively. 
Purple, pink and green lines represent the shear, F flexion and G flexion 
signal if the halo is a smooth SIS (dashed). The solid lines are an SIS pro- 
file as fitted to the shear data points in a simulated galaxy containing TSIS 
subhaloes and translated into predictions for the flexions. 



5.5 The Effect of Substructure 

Since flexion is more sensitive to the underlying mass distribution 
on small scales than shear is, we expect it to respond differently 
to the presence of substructure in galaxy haloes. To test whether 
this has any impact on our analysis we take a galaxy- size SIS halo 
(see B06, for shear and flexion expressions) and populate it with 
subhaloes, allowing 20% of the mass to be in substructure. The 
total mass of the halo is 10^^ h~^MQ and the galaxy is placed 
at ^ = 0.35 with Di/Dis = 0.5. We spread the substructure 
mass over 100 subhaloes, randomly distributed according to an SIS 
density profile. Finally we average the azimuthally averaged sig- 
nal over 100 such galaxies. Now, subhaloes are generally stripped. 
To approximate this we use a Truncated SIS (TSIS) profile for the 
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Figure 7. The correlation matrix between the shear and flexion bins, using 
5000 bootstrap resamples. Please note the scale; to display the minute vari- 
ations between off-diagonal elements we have artificially set the diagonal 
elements (dark green) only to 0.1, whilst all other elements are unsealed and 
normalised to diagonal elements of 1.0 as is customary. 

subhaloes fsee lHoekstra et alJl2004 for constraints on parameters). 
The TSIS convergence is given by 

^(0) = ^ f 1 - -r^=\ (17) 

where Oe 'i^ the Einstein radius and Os is a truncation scale where 
the profile steepens. On small scales {0 <^ Os) the TSIS behaves 
like an SIS but at large sc ales (9 ^ Os) the profile decreases as 0^. 
The TSIS shear is given in l Schneider & Rixl (Il997h and the flexions 
are 

^^'^ = w{wTW^-')^'' '''' 

and 

where is the position angle of the background source. Using the 
parameters above and a truncation scale = 2 arcsec for the sub- 
haloes we get the results shown in Figure [6l The shear profile fit is 
pulled down slightly compared to a smooth halo but the flexions 
are not similarly affected. Due to the substructure the flexions are 
more scattered, but the overall trend is for the points to follow the 
smooth profile, or even slightly above in the F flexion case. Thus 
the flexions seem overestimated compared to the shear fit. We stress 
however that the fraction of substructure used in this test (20%) is 
high to exaggerate the effect. The test does show that substructure 
may affect the flexions differently to the shear, but its influence is 
likely less than the excess currently observed in COSMOS. 

5.6 Profile Determination 

One of the most interesting potential uses of flexion is as an aid 
to shear in determining the inner density profiles of dark matter 
haloes. The two signals are sensitive to the underlying density pro- 
file on different scales, so combining the two will give us tighter 
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Figure 8. Joint profile constraints using shear and F flexion. The top (bot- 
tom) panel shows the results for the power-law (NFW) fit. Purple (thin solid) 
lines represent shear and green (dashed) represent F flexion. The contours 
show the 67.8%, 95.4% and 99.7% confidence limits respectively in terms 
of constant Ax^ (2.30, 6.17 and 11.8 respectively). The white (thick) con- 
tour marks the joint confidence limits. The grey-scale is logarithmic in x^- 



constraints than either on their own. To combine the shear and flex- 
ion signals we have to take any correlation between them into ac- 
count. B06 assumed that the shear and flexion measurements would 
be uncorrected. Here we confirm this assumption through the cor- 
relation matrix between the shear and flexion bins, using 5000 boot- 
straps, shown in Figure [T] This implies that it is trivial to combine 
the shear and flexion information to find the profile of an average 
lens. We use the F flexion in conjunction with the shear to fit density 
profiles to the measured signal. For this purpose we try two differ- 
ent families of profiles: the power-law and the NFW. Our general 
power-law is defined as 

7 - -Ad-"" (20) 

with d the distance from the lens, and the amplitude A and the 
index n free parameters. An index of n = 1 would be equivalent to 
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an SIS. The above expression is easily differentiated to give the F 
flexion 

{n-2)Ad-''-\ (21) 

The expressions for the NFW profi les are somewhat more com pli- 
cate d but they are given in full in Wright & BrainerdI (1200 0') and 
Ib06 for shear and flexion respectively. Here we leave the virial ra- 
dius M200 and the concentration c as free and independent param- 
eters. We fit the power-law and NFW profiles to the inner 100 kpc 
only as this is the region where F flexion becomes important and 
the shear profile is not affected by halo-halo contamination. 

The top panel in Figure [8] shows that both the shear and the 
F flexion are consistent with an SIS (n = 1), although together they 
prefer a slightly lower power-law index of n = 0.731q 43. The bot- 
tom panel shows that it is difficult to constrain the NFW concentra- 
tion if it is left completely unrestricted. This analysis with two free 
and independent parameters is not completely representative, how- 
ever, since s imulations indicat e a fixed mean mass-concentration 
relationship (iDuffv et ahkOO SV It is also important to note that the 
average profile we constrain here is a composite of lenses in a large 
redshift range. Detection at the high end of the redshift distribution 
tend to be biased towards intrinsically brighter objects than at the 
low end. We also combine measurements from lenses of different 
sizes and morphologies. Nonetheless, combining shear and F flex- 
ion does provide tighter constraints than shear alone on the density 
profiles, and this is an important proof of concept. The resulting 
mass estimate for the average lens in COSMOS from the combined 
NFW fit is M200 = 2.12l5;g° x 10^^ H-^Mq with a concentra- 
tion of c = 4.82t^;^. 



6 DISCUSSION AND CONCLUSIONS 

We have shown a significant detection of galaxy-galaxy F flexion 
for the first time with Shapelets using the space-based COSMOS 
data set. We used this flexion signal in conjunction with the shear 
to constrain the average density profile of the galaxy haloes in our 
lens sample. We found a power-law profile consistent with an SIS. 
Furthermore, we showed that the inclusion of F flexion provides 
tighter constraints on both power-law and NFW profiles, an impor- 
tant proof of concept. 

The galaxy-galaxy F flexion signal measured in COSMOS is 
slightly higher than expected from the shear signal, especially if we 
apply the multiplicative bias correction. There is however no indi- 
cation from the cross term that there are systematics present. The 
discrepancy could be partly due to insufficient nearby object light 
removal, but this is unlikely to explain the full offset. Substructure 
in galaxy haloes may cause excess F flexion compared to what the 
shear measures. However, a large fraction of the galaxy halo mass 
has to be in substr ucture in order for the effe ct to become signif- 
icant. We note that iGoldberg & BaconI (l2005h also find shear and 
F flexion signals that are inconsistent with each other; the velocity 
resulting from an SIS profile fit to their F flexion signal is nearly 
twice that found using shear. This is qualitatively consistent with 
our findings, which leads us to believe that there is something more 
fundamental affecting the signal. In the near future we would like 
to further investigate the dependence of these discrepancies on lens 
properties. 

We measure a galaxy-galaxy G flexion signal that is consistent 
with the predicted profile, but due to the large measurement errors 
it is also consistent with zero. This measurement is a lot noisier 
than the other two, an effect most likely caused by the fact that 



there is less information available in the higher m-order Shapelets 
for fainter sources. To measure a G flexion signal we thus require 
many well-resolved sources, an extravagance not yet awarded us. 
Future large space-based surveys such as EUCLID will enable us 
to investigate G flexion further, but for now F flexion is a promising 
tool in its own right. 

The software introduced in this paper, the MV pipeline, is able 
to detect these higher order lensing distortions. We have shown that 
in practice, the Shapelets F flexion measure is affected by light from 
nearby bright objects and detailed a way to correct for this effect. 
This BOR does require further sophistication to account for large, 
well resolved galaxies, galaxies which are not well described by 
the single Sersic light profile employed here. From the FLASHES 
simulations it is clear that there is more work required in order to 
improve the accuracy of the F flexion measurements for future sur- 
veys. Noise related biases are particularly significant for this type 
of shape measure, and we have modeled these biases in COSMOS. 

In the future we hope to measure flexion on a larger survey, 
enabling us to reduce the noise so that we can investigate the trend 
with e.g. redshift and lens mass. A larger number of sources would 
also enable us to further tighten the profile constraints in the inner 
regions of dark matter haloes where baryons become important. It 
is not yet clear how well we can measure flexion on ground-based 
data, but surveys like KiDS, CFHTLS and RCS2 should provide an 
excellent test-bed. 
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Figure Al. The multiplicative bias on the second component for each of 
shear, F flexion and G flexion. The purple stars represent shear, pink circles 
represent F flexion and green triangles represent G flexion. The symbols 
and solid lines show the weighted averages whilst the dashed lines show the 
CHP average. 
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Figure A2. The additive bias on the first component for each of shear, F 
flexion and G flexion. As before, the purple stars represent shear, pink cir- 
cles represent F flexion and green triangles represent G flexion. The sym- 
bols and solid lines show the weighted averages whilst the dashed lines 
show the CHP average. 



APPENDIX A: FLASHES RESULTS 

The figures shown in this Appendix complement Figure [2 in the 
main paper and provide additional detail on the results from run- 
ning the MV pipeline on FLASHES, with rimax = 10. The param- 
eters m and c are defined through 

/ measured\ input input , /ai\ 

and similarly for the flexions, where i = 1, 2 is the component. We 
use two different techniques to estimate the average distortion on 
each image: a weighted average and Convex Hull Peeling (CHP). 

In Figure lAll we show the multiplicative bias of the second 
component for each of shear, F flexion and G flexion as a func- 
tion of the different simulation branches. For these results we use 
a Shapelets order of rimax = 10. This bias behaves as the multi- 
plicative bias of the first component (Figure |2]), as expected. The 
biases of all distortion measurements, and in particular F flexion, 
are severely dependent on S/N and brightness profile. 

The additive bias c is minimal for shear and F flexion (see 
Figures lA2l and lA3l l indicating that the PSF is either well corrected 
for or not significantly affecting these two measurements. For the 
G Flexion the offset is larger. 



APPENDIX B: COSMOS DATA ANALYSIS 
Bl Catalogue Creation 

To maximise the number of lens- source pairs we use all objects 
with assigned photometric redshifts as sources, but imposing a red- 
shift cut of z < 0.6 for lenses. Additionally we use sources without 
individual redshifts (SIO redshift bin 6), assigning mean angular di- 
ameter distance ratios (Ds/Dis) to these lens-source pairs accord- 
ing to the estimated redshift distribution employed by SIO. We then 
weight all pairs with their individual lensing efficiency, similar to 



the weighting scheme in e.g. iMandelbaum et al.l (l2006bh (see Ap- 
pendix |B4l). This downweights pairs that are close in redshift and 
naturally removes pairs where the "source" is in front of the "lens". 
To the source catalogues we apply the following cuts: 

• S/N > 10. This cut is important as the F flexion measure- 
ment in particular gets heavily biased towards low S/N (see Sec- 
tion l432l) . 

• If the centroid cannot be determined accurately the Shapelets 
decomposition will be inferior. Therefore objects where the code 
is forced to move the centroid compared to the one estimated by 
SExtractor by more than half a pixel are excluded. 

• The summed power in constant m of the polar Shapelets pro- 
vides an indicator of the Shapelet fit being affected by a neighbour- 
ing object. If the fractional power is particularly high at high orders 
the object is excluded (see KK06, for more details). 

• If the FWHM or scale radius of the object is too small com- 
pared to the scale radius of the PSF the object is excluded. 

• If 7^ > 1.4, J^^ > 3.0 arcsec"^ or > 6.6 arcsec"^ then 
the object is excluded. These numbers are based on the measured 
distributions and the cuts are applied to remove outliers with very 
noisy shape measurements. 

• Finally, we remove faint objects with an assigned photometric 
redshift of ^ < 0.6 that h ave a prominent secondary peak at Z2nd > 
0.6, as discussed in lSlOl 

B2 PSF Interpolation 

The ACS PSF fluctuates both spat i ally and temporally (e.g. 
[Rhodes et aD l2007l : ISchrabback et al.l l2007h . a variation mostly 
driven by changes in telescope focus caused for example by the 
breathing of the telescope. We can map the PSF using stars, but in 
high-galactic latitude ACS fields typically only ^ 10 — 20 stars 
are present. This number is too low for the standard approach of a 
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Figure A3. The additive bias on the second component for each of shear, F 
flexion and G flexion. As before, the purple stars represent shear, pink cir- 
cles represent F flexion and green triangles represent G flexion. The sym- 
bols and solid lines show the weighted averages whilst the dashed lines 
show the CHP average. 



poly nomial interpolation. Instead, we closely follow the analysis of 
ISIQI who conducted a principal component analysis (PCA) of the 
ACS PSF variation as measured in dense stellar fields. They found 
that ^ 97% of the PSF variation can be described with a single 
parameter (the first principal component). This parameter is related 
to the HST focus position, and we therefore dub it 'focus {3 

Here we make use of the SIO measurement of the HST focus 
in all COSMOS exposures and the investigated stellar field expo- 
sures. We also obtain Shapelets versions of the focus-dependent 
lSIO PSF models, by decomposing the dense stellar field stars into 
Shapelets and interpolating between them with polynomials which 
are varied both spatially and with different powers of the focus 
principal component coefficient. From these models and from the 
COSMOS focus estimates we then compute a Shapelets PSF model 
for each COSMOS exposure, which we then combine to obtain a 
model for the stacked PSF at all galaxy positions. 



B3 CTI Correction 

Our correction for CTI again closely follows SIO, who derive para- 
metric corrections for the change in polarization for both galaxies 
and stars. The correction for stars is important in order to measure 
the actual PSF, independent of the non-linear CTI effects. In the 
stellar field analysis we therefore correct the PSF cartesian Shapelet 
coefficients for CTI before generating the PCA PSF model. In or- 
der to estimate the influence of CTI on the different Shapelet coef- 
ficients, we follow SIO and spatially fit each coefficient within one 
exposure with polynomials. Due to the limited depth of the charge 



^ The capturing of small additional variations beyond focus was relevant 
for the cosmic shear analysis of SIO. Here we can safely ignore these mi- 
nor additional effects. Galaxy-galaxy lensing is much less sensitive to PSF 
anisotropy residuals as they cancel out to first order. 
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Figure Bl. CTI-induced residuals on the stellar Shapelet coefficients S20 
(left) and so3 (right) in an example star field. The black stars show the 
mean of the coefficients as a function of stellar flux after subtraction of a 
spatial third-order polynomial model derived from bright stars to separate 
PSF and CTI effects. Each coefficient has been scaled to a reference number 
of ytrans = 2048 parallel readout transfers. The purple curves show the 
parametric CTI model, jointly determined from 700 stellar field exposures. 
The horizontal dashed line indicates an offset corresponding to the mean 
CTI model for the bright stars used for the polynomial interpolation. The 
green triangles indicate the corrected coefficients after subtraction of the 
CTI model. 



traps, CTI is non-linear, and has a larger relative impact on faint 
sources than on bright ones. The CTI effect can thus be estimated 
from the flux-dependent residuals, after the polynomial model has 
been used to subtract both the flux-independent PSF and the flux- 
averaged CTI signal. 

Figure IB II shows these residuals as a function of stellar flux 
for the stellar Shapelets coefficients S20 and sos in one example 
stellar field. Here the residuals were scaled to the same number 
of readout transfers (2048). The CTI effect on the coefficients is 
clearly visible (black stars), but with our power law model (curve) 
it can be well corrected for (green triangles). The model is fit si- 
multaneously from all 700 stellar fields as a function of stellar flux, 
sky background, time and number of readout transfers (see SIO). 
CTI affects object shapes in the readout direction, which also af- 
ter drizzling roughly matches the ^/-direction. Thus CTI residuals 
are expected to be roughly symmetric about the y-SLxis and hence 
vanish for coefficients Sab with odd a. In the drizzled images the 
readout direction is up for the upper and down for the lower chip 
and the CTI trails occur in the opposite directions. This leads to a 
sign switch for coefficients with basis functions that are not sym- 
metric about the x-axis (odd 6), and we have taken this into account 
for so3 in Figure IbT] We have detected (and modeled) a significant 
signature of CTI on the following stellar Shapelets coefficients: soo, 
S02, so3, so4, so5, S20, S21, S22, S40, and seo- 

The correction of galaxy shapes for CTI again closely follows 
SIO. Here we fit power-law corrections to the shear and (now in 
addition) flexion estimates as a function of galaxy flux, flux radius, 
sky background, tim e, and number of readout- transfers. Note that 
iMassev et all (l2010h introduced a more advanced CTI correction 
scheme operating directly on the pixel level. This is expected to 
yield higher precision, enabling for example the correction of the 
soi component, which cannot be estimated with our method due to 
its degeneracy with a simple shift in object position. However, we 
are confident that our correction scheme is sufficiently accurate for 
the analysis presented here, in particular as potential residuals can- 
cel to first order for the azimuthally averaged galaxy-galaxy lensing 
signal. 



© 2010 RAS, MNRAS 000,[T}{T3] 



Flexion in COSMOS 13 



B4 Signal Computation 

For the Navarro - Frenk - White (NFW) 
(iNavarro. Frenk. & Whitd Il996h . the strength of 
signal scales as 
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where Du Ds and Dis are the angular diameter distances to 
the lens, to th e source, and between lens and source respectively 
dWright & Br ainerd 2000). The flexion signals scale as 
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We therefore weight the signals accordingly, scale them to a 
reference lens and source redshift and compute the weighted aver- 
age in 25 logarithmic distance bins as follows: 
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and similar for the flexions, with the shear estimator and weight 
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is the geometric lensing efficiency and a^^i = cr^ 
c^^2,^ is the error on the shape measurement with cr^, intr the intrin- 
sic shear noise. By contrast we use the following F flexion estimator 
and weight: 
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and similarly for the G flexion. 
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APPENDIX C: HIGH REDSHIFT RESULTS 

As specified in the main paper, the lens catalogue we use has a 
redshift cut of z < 0.6. This is to avoid having to go too close to the 
lens on the sky in order to see a flexion signal. Within an angular 
radius of 2 arcsec we have low confidence in the results; we are 
simply too close to the lensing galaxies and it becomes difficult to 
account for effects induced by the lens light. BOR corrects for light 
leakage at larger radii, but the correction is most likely incomplete 
very close to the lens due to deviations from a smooth Sersic profile. 
For objects beyond our lens sample, the median redshift is close 
to 1.0. At this redshift the angular distance limit of 2 arcsec on 
the sky corresponds to a physical distance of about 17 kpc. The F 
flexion falls off to low values already at about 20 kpc for a typical 
galaxy, so we are left with a very low signal within a narrow ring 
around the lens. Imposing the redshift cut of z < 0.6 on lenses 
gives us a median lens redshift of z = 0.27 at which the inner limit 
corresponds to 9 kpc, leaving a wider distance interval in which 
we can investigate the F flexion signal. 

In Figure ICll we show the galaxy-galaxy signal for the high 
redshift sample, i.e. for lenses with z > 0.6. The bins that are 
within 2 arcsec of the average lens in this sample, and which are 
most likely contaminated by lens light, are marked with dotted 
lines. The F flexion signal outside of this limit does agree well with 
the profile predicted by the shear, but falls off quickly. 
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Figure Dl. A comparison between the shears used in this paper and the 
ones used in SIO. Black stars (green circles) show the difference between 
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APPENDIX D: COMPARISON WITH KSB 

We compare our galaxy-galaxy shear signal to the one we get using 
the shears from . S 10 , using all the cuts normally applied in each 
analysis so that only common objects are used. The bias correction 
described in their paper is incorporated in their shears whilst our 
measurements have no correction applied. However, due to our S/N 
cut (see Appendix IB 11) their correction is always less than 4.2%. 

As shown in Figure lPTI the difference between the results from 
the two pipelines, KSB and Shapelets, is negligible. This provides 
an independent confirmation that the MV pipeline produces shears 
of as high a quality as the state-of-the-art weak lensing analysis 
presented in lSld . 



© 2010 RAS, MNRAS 000,[T}tT3] 



14 M. Velander, K. Kuijken & T. Schrabback 




10 20 50 100 200 500 10 20 50 100 200 500 10 20 50 100 200 500 

Distance [kpc] Distance [kpc] Distance [kpc] 



Figure CI. The galaxy-galaxy lensing results from running the MV pipeline on the COSMOS data, with Timax — 10. Black solid points represent the 
tangential signal and green circles represent the cross term. Empty circles with dotted error bars are bins that are too close to the lens on the sky. Please note 
that the SIS and NFW profiles have been fitted to the shear data and then translated into predictions for T and Q curves. 
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